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Abstract 

The CP conserving amplitudes for the decays K — * 3n are calculated in Chiral 
Perturbation Theory. The calculation is made at the next-to-leading order including 
strong and local electromagnetic isospin breaking. A comparison is made between 
the squared amplitudes with and without isospin violation to estimate the size of the 
effect. We find corrections of order five percent in the amplitudes. 



PACS numbers: 13.20.Eb; 12.39.Fe; 14.40.Aq; 11.30.Rd 



1 Introduction 



Chiral Perturbation Theory (ChPT) is the effective field theory of the low-energy strong 
interactions. It was introduced in its present form by Weinberg, Gasser and Leutwyler 
[Tj 12 El and it has had many successes and applications. A pedagogical introduction can 
be found in [4j. The theory has been applied as well for nonleptonic weak decays. The 
main work of extending ChPT to the nonleptonic weak interaction was done by Kambor, 
Missimer and Wyler who worked out the general formalism and applied it to K — > 27r, 37r 
decays [SJ. These results were then used to obtain directly relations between physical 
observables in jjj. Reviews of applications of ChPT to nonleptonic weak interactions are 
jH]. The expressions of Ref. [S] were never published and have been lost. 

The next-to-leading order K —>■ 2n amplitudes in the isospin limit were recalculated in 
pi ITU] and the K — > 3tt amplitudes in ^Oj- The latter have been confirmed in [TJ]. In [10] 
a full fit to all available K — » 2tt and K — > 3ir experimental results was performed and it 
was found that one could fit the decay rates and the slopes in the Dalitz plot. There was 
a discrepancy with the observed quadratic slopes or curvatures in the Dalitz plot and this 
can have several different origins. It could be an experimental problem (especially given 
the discrepancies between several experiments) or it could have a theoretical origin. In 
the latter the amplitudes calculated so far have three different types of corrections, isospin 
breaking, electromagnetic corrections or higher order ChPT corrections. In this paper we 
investigate the first class. Work is in progress to evaluate the electromagnetic corrections 
as well. More precisely, in this paper we investigate effects of the quark mass difference 
m u — r^d as well as the local electromagnetic effects. We have worked out all orders 
in m u — rrid to order four in the chiral expansion but we have found that the numerical 
difference between the all order isospin breaking and the first order effect was much smaller 
than the other uncertainties in the calculations. We therefore present results only to first 
order in m u — m<j. 

There exists very little work on isospin breaking in K — > 3ir decays and none that we 
know of within ChPT. The K — > 2tt case has been investigated more extensively since it 
has possibly strong effects on e'/e [TJ]- Recent works are [13] ITU IT5| ITS|. 

The outline of this paper is as follows. The next section describes strong isospin break- 
ing. Section[3]presents the Chiral Lagrangians needed for the calculation to next-to-leading 
order and section 0] specifies the decays and describes the relevant kinematics. The analyt- 
ical results, the amplitudes for K — > 3n, are described in section [5] and in section [0 some 
numerical results are shown. The last section contains the conclusions. 

2 Strong Isospin Breaking 

Historically, the isospin approximation meant treating the proton and neutron in the same 
way. In a quark picture this corresponds to treating the up- and down-quark as being 
identical. This implies putting m u = rrid as well as neglecting electromagnetism. Obviously 
this is an approximation but since one can roughly estimate the corrections, and it greatly 
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simplifies the calculations, its use is very common. However, to reduce the errors one has 
to quote, the isospin approximation has to be abandoned. 

Isospin breaking is usually divided into two parts. Strong isospin breaking coming from 
the fact that m u ^ ma, and electromagnetic isospin breaking from including electromag- 
netic effects. In this paper we only take into account the strong part and the local part 
(the part which doesn't include explicit photons) of the EM isospin breaking. From now 
on we will use the name strong isospin breaking for the sum of these two effects. 

The difference between m u and leads to mixing between tt and rj. This means 
changes in the formulas for both the physical masses of tt and rj as well as the amplitude 
for any process involving either of the two. For a detailed discussion see [T7j . 

Including the local EM part means introducing new Lagrangians at each order, propor- 
tional to e 2 and p 2 e 2 respectively. 



3 The ChPT Lagrangians 

The starting point of our ChPT calculation is the Lagrangians. The order parameters in the 
perturbation series are p, the momenta of, and m, the mass of the pseudoscalars. Including 
isospin breaking also e, the electron charge, is used as an order parameter. Leading order 
(order two) therefore means terms of order p 2 ,m 2 ,e 2 , and next-to-leading order (order 
four) p A , p 2 m 2 , m 4 , p 2 e 2 and m 2 e 2 (e 4 is neglected). 



3.1 Leading Order 

Our leading order effective Lagrangian assumes the form 

£2 = £>S2 + £>W2 + £>E2 



(1) 



Here £52 refers to the strong AS = part, C W2 the weak AS = ±1 part, and Ce2 the 
strong-electromagnetic and weak-electromagnetic parts combined. For the strong part we 
have 

C S 2 = ^ <«X + X+> (2) 

Here (A) stands for the flavour trace of the matrix A, and Fq is the pion decay constant 
in the chiral limit. We also define the matrices u^, u and x± 



■ ~ - U , x± = u ] xu ] ± ux ] u . 
where the special unitary matrix U contains the Goldstone boson fields 



(3) 



U = exp 



M 



( >3 + > K+ \ 



(4) 
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The formalism we use is the external field method of 0, but for our purpose it suffices to 
set 



V 



m d 



and D„U = d^U. 



(5) 



We diagonalize the quadratic terms in Q by a rotation 

7T° = 7r 3 cos e + 7/s sin e 
rj = — 7r 3 sin e + r/ 8 cos e , 

where the lowest order mixing angle e satisfies 

y/3m d - m u 



tan(2e) 



2 m,-m 



(6) 



(7) 



with 2m = m u + m^. 

For the weak part the Lagrangian has the form ^Hj 



C 



W2 



+ h.c. 



The tensor t %3> has as nonzero components 

1 



.21,13 _ +13,21 _ ± . +22,23 _ +23,22 

3 ' 
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23,33 _ +33,23 



t 



and the matrix A,-,- is defined as 



— ; t 
6 ' 



23,11 _ ^11,23 



Ajj = uXijW , (Ay) o6 = 5j& . 

The coefficient C is defined such that in the chiral and large N c limits G$ = G?i = 1, 



(9) 



(10) 



11^ 



Finally, the remaining electromagnetic part, relevant for this calculation, looks like (see 
e.g. HH) 

C E2 = e 2 F*Z(Q L Q R ) + e 2 F*(TQ R ) (12) 
where the strangeness-conserving term includes 



Ql = uQu) , Q R = u^Qu with Q 



2/3 
-1/3 
-1/3 



(13) 



and the weak-electromagnetic term is characterized by a constant Ge (Ge = gewkGg in [T9]). 

T = G s F 2 A 32 + h.c. (14) 
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3.2 Next-to-leading Order 

Since ChPT is a non-renormalizable theory, new terms have to be added at each order to 
compensate for the loop-divergences. This means that the Lagrangians increase in size for 
every new order and the number of free parameters rise as well. At next-to-leading order 
we split the Lagrangian in four parts which, in obvious notation, are 



£ A = £ 



S4 



£ 



Wi 



£ 



S2E2 



+ £ 



W2E2 



(G 8 ). 



(15) 



Here the notation (G 8 ) indicates that here only the dominant Gg-part is included in the 
Lagrangian, and therefore in the calculation. The strong part we need looks like [21 [20] 

£s4 = L x (u/) 2 + L 2 (it/) + L 3 (u^u u u v ) + U («/) (x+> 

+ L 5 + U ( X+ ) 2 + L 7 (x-) 2 + L 8 i ((xl) + (X-) ) , (16) 

and the weak part, quoting only the terms relevant for K — > 3tt decays, is l2*T| I22j. 
£wa 



C Fq < Gf. 



NiOf + N 2 0\ + N 3 8 3 + NJDl + N 5 8 + N 6 8 + N 7 O s 7 



+N 8 8 8 + N 9 O s 9 + N 10 O 8 l0 + N u 8 n + N l2 8 2 + N 13 



18 



+G- 



27 



D t Of + D 2 Of + D 4 Of + D 5 Of + D 6 Of + D 7 Of + D 2e O 



\27 



\27 



27 



\27 



127 



127 

26 



+DvtO% + D 28 Of s + D 29 Of + D 30 OH + D 31 0*l 
The octet operators are 



+ h.c. . 



(17) 



o 8 = 


(A 32 u^u u u u ) 


o 8 = 


{A 32 u^u u u v u^) 


o 8 3 = 


(A 32 u^u u )(u»u u ) 


o 8 = 


(A 32Ufl )(u u u»u u ) 


oi = 


(A 32 (x+V^ + w/x+)) 


oi = 


(A 32 n M )(^ X+ ) 


o 8 = 




o 8 - 

^8 — 




o 8 - 


(A 32 (x-w/-it/x-)) 


o 8 - 

^10 ~ 


; (A 32 x+X+) 


o 8 - 


: (A 32 x + )(x + ) 


o 8 - 

w 12 - 


: (A32X-X-) 


o 8 - 

^13 ~ 


= (A 32 x-)(x-), 







and the 27 operators are 

Of 
Of 

of 
of 



t^ kl (A ijX+ )(A klX+ ) 
t ij > kl (A ijX -)(A klX -) 
t ij > kl (A ijU „)(A kl (u>*x + + X+un) 
t^ kl (A l3 u,)(A kl (u»X--X-un) 
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of 


= t ijM 


(A ij x+)(^kiUuU fl ) 


of 


= t ij ' k 


(A.ijUn) (Afc/w^) (x+) 


O 27 

w 26 


— tf>3,k 


{A lJ u^){A kl u v u u ) 


Of 


= j-ij,k 


(Aij(u^u u + u v u tl )){A k iu tl u v ) 


of 


= f j ' k 


(A i:j (u^u u - u u u ll )){A k iu^u v ) 


O 27 

^29 


= f j ' k 


(A lj u fJ )(A kl u u u> M u l/ ) 


O 27 




(A ij u ll )(A kl (u'*u v u v + m/u")) 


O 27 


= t ij ' k 


{A l3 u^){A kl u^){u u u v ) . 



The complete minimal Lagrangian of 0(G&e 2 p 2 ) takes the form ^H] 

C W2E2 (G 8 ) = G 8 e 2 F*J2 z iQt + h.c. 



(19) 



(20) 



i=l 



with operators Q™ of 0(p 2 ) and dimensionless coupling constants Zi. A linear independent 
set of operators is given by 



Qr = 


(A 32 {Q R , X+ }) 


Q w 2 = 


(A 32 Q R )( X+ ) 


Ql" = 


(A 32 Qr)(x+Qr) 


Qt = 


(A 32X+ )(QlQr) 


Q% = 


(A 32Ufl u») 


Q w ,= 


(A 32 {Q R ,u^}) 


Q w ,= 


(A 32 u^)(Q L Q R ) 


Qs = 


(A 32 u^)(Q L u») 


Ql" = 


(A 32 u^(Q R u») 


Qio = 


(A 32 ^)({Q L ,Q /? K) 


Qfi = 


(A 32 {Q R ,u^})(Q L u») 


Ql2 = 


(A 32 {Q R , Utl })(Q R u») 


Ql3 = 


(A 32 Q,)(n/) 


Ql4 = 


(A 32 Q R )(u^Q R ) . 



(21) 



Finally the strong-electromagnetic part, which looks like [2lH 12*3] 

n 

Cs2E2 = e 2 F 2 J2K i Q s l 



(22) 



i=l 



with 



Qs = ((QluJ (Qlu») + (Q R ^) (Q r vT)) 

Ql = ({Ql + Ql)u^) 
Q s 7 = \(Ql + Ql) <x+> 

Ql = ({Qt + Ql) x+ ) 

Qli = ((QrQl — QlQr)x~) ■ 



Q s 2 
Ql 
Ql 
Ql 
Q\o 



(QlQr) 
(Ql<v) (QrU») 

((QlQr + QrQi>X> 
(QlQr) (x+> 
= ((QlQr+QrQl)x+) 



(23) 
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7. 




K 




i 
i 


1/ / Z4 + 9/ Z Z — 9/4 Pewk 


i 
i 


3 /S 
— d / O 


o 
L 


1/9 8/^7 1 /On , 

— i/ z — e/o Ai — i/ z p ew k 


9 
Z 


1 /9 7 

— 1 / Z Z/ 


Q 
O 


q /c 7/9 7 

— O / O — f / Z Z/ 


Q 
O 


o / o 


A 


1 /8 _J_ 7 /9 7 

9/ O "T 1 / Z Zj 


A 


7 
— Z 





1 
1 


o 


Q /S 
y / o 


u 


7/4 ^ /9 7 — 3 /A n , 


u 


^ /4 7 

— O / ft Z 


7 




7 


n 


8 


1/4 


8 


-1/2Z 


9 


11/12- 2/3 Z-^ewk 


9 


1/8 


10 


3/4+ 1/2 zT 


10 


-1/8- 3/4 Z 


11 


3/4 + zT 


11 


-1/16 


12 


-3/4 






13 


35/24 + 3/2 Z- 1/2 ^ ewk 






14 


-3/2 - 15/2 zT 







Table 1: The coefficients of the subtraction of the infinite parts defined in Eq. ([24)1 . 

The infinities appearing in the loop diagrams are canceled by replacing the coefficients 
in ([16)1 . ()17|) . ()20|) and (|2*2*|) by the renormalized coefficients and a subtraction part. The 
infinities needed in the strong sector were calculated first in [2j and those for the weak 
sector in [Sj. The terms are all of the type 

X t = (e^ (XI + . (24) 

with the dimension of space-time d = 4 — 2e and 

c=-i(ln(47r)+r'(l) + l) ■ (25) 

For X = L, N, D see [H], and for X = Z, K see Table □ O EBj- However, since the full 
e 2 p 2 -contribution is not included (the explicit photon diagrams are left out), the infinities 
don't cancel completely. They still provide a number of checks on the calculations though. 
Since contributions of order e 4 are neglected, Lagrangians allready containing e 2 will not 
add any explicit photon terms. Therefore the infinities proportional to Z and Ge cancel. 
Even more, since there are no explicit photon contributions to the decays involving only 
neutral particles, like Ki — > 7r 7r 7r°, the infinities in that channel cancel completely. 

4 Kinematics 

There are five different CP-conserving decays of the type K — > 37r: 

K L (k) - APi)AP2)APz), [^ooo], 
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K L {k) -> 7r + ( Pl )7r-(p 2 )7r (p 3 ), [^_ ] , 

K s (k) -> 7r + ( Pl )vr-(p 2 )7r (p 3 ), [Aj_ ] , 

#+(*) _> ^^^(pajTr+Cps), [A 00+ ], 

- 7r + ( Pl )vr + (p 2 )7r-(p3), , (26) 

where we have indicated the four-momentum defined for each particle and the symbol used 
for the amplitude. The K~ decays are not treated separately since they are counterparts 
to the K + decays. 

The kinematics is treated using 

s 1 = (k-pif, s 2 = (k-p 2 ) 2 , s 3 = (k-p 3 ) 2 . (27) 

The amplitudes are expanded in terms of the Dalitz plot variables x and y defined as 



s 3 ~ S S 2 - Si 1 

V = o — ' x = 2 — ' s = -{s 1 + s 2 + s 3 ) . {21 



The amplitudes for Kl — > 7r + 7r~7r°, K + — > tt + 7t + tt~ and K + — > 7r°7r°7r + are symmetric 
under the interchange of the first two pions because of CP or Bose-symmetry. The ampli- 
tude for Ki — > 7r°7r 7r° is of course symmetric under the interchange of all three final state 
particles and the one for K$ — > 7r + 7r~7r° is antisymmetric under the interchange of 7r + and 
7r~ because of CP. 

The amplitudes to order four can be written in terms of single variable functions Mj(s) 
[E] as 

A^ ( Sl ,s 2 ,s 3 ) = Afo(si) + M (s 2 ) + M (s 3 ) , 

Aj_ (si, 8 2> s 3 ) = M 1 (s 3 ) + M 2 (si) + M 2 (s 2 ) + M 3 (si)(s 2 - s 3 ) + M 3 (s 2 )( Sl - s 3 ) , 
Aj_o(si,S2,s 3 ) = M 4 ( Sl )-M 4 (s 2 ) + M 5 ( Sl )( S2 - S3 )-M 5 (s 2 )( Sl -s 3 ) 
+M 6 (s 3 )(si-s 2 ), 

s 2 , s 3 ) = M 7 (s 3 ) + M 8 (si) + M 8 (s 2 ) + M 9 ( Sl ){s 2 - s 3 ) + M 9 (s 2 ){ Sl - s 3 ) , 
s 2 , 8s) = M 10 (s 3 ) + Mii(si) + M u (s 2 ) + M 12 ( Sl )(s 2 - s 3 ) + M 12 (s 2 )( Sl - s 3 ) . 

(29) 

The functions Mj(s) are not unique. The full list of ambiguities can be found in App. A 

of mg. 



5 Analytical Results 

5.1 Lowest order 

The two diagrams that contribute to lowest order can be seen in Fig^ 

One of the problems in expressing the amplitudes is the various possible choices of 
masses to express the isospin conserving piece in. The simplest expressions usually follow 
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Figure 1: The tree level diagrams for K — > 3n. A filled square is a vertex from Cw2 or 
Ce2{AS = 1) and a filled circle from Cs2 or Ce2{AS = 0). 

when using a mixture of the neutral and charged masses and in the part that is obviously 
isospin violating we can use the neutral and charged masses interchangeably, denoted by 
m\ and ml- 

The simplest expressions for the lowest order we have found are 

M (s)\ p , = i (G 8 - G 27 ) \m\, (l + y/3 sine) . (30) 

\ i . CFn f , „ _ v 1 o sine 13 \ 9 

Mi( S )| p2 = tj^UG 8 -G 27 )-m 2 K0 + -^[G 8 --G 27 jm 2 K 



+Ze 2 F 2 -G 8 + 



_G | G 27 6rn? K - 16ml ' 
3 8 3 mi- — ml 



M 2 (s)\ n 2 = 0, 



CFlj_Gs (11 ml + 1 mlo \ 

F*F K \ 3 + ^ 27 { 2 + 3m 2 K+ -ml+ 2m 2 K0 -ml ) 

F— 7=- f-2G 8 + + fz Gs + 1^ -^U. (31) 



\/3 \ 6 - ml J \ 2 



M 4 (s)| p2 
M 5 (s)| p2 

M 6 ( S )| p2 = z^l^|G 27 

-G 8 - 




5 

2 ~ 

G 2 7 7m^ 




This is in the isospin limit the same as (31) in JU] using Eq. (A. 9) there. 



S 



Mgfs)!^ = 0, 

+ — ( Gs + %I4^U(Z G8+ 1 GE )^U. (33) 



\/3 \ 3 mff - ml J \ 2 J m 2 K - ml 

This can be brought into the form of (32) of ^0] using the equivalents of (A.3-A.5) there. 

/ \ i ■ [ „ / 2 o \ / 13 n 13 2 

Mio(s)| p 2 = 2-^|G 8 - m x+ - m T+ J + G 27 (^-ys + m x+ + ym r 

-8Ze 2 F 2 G 8 - fze 2 F 2 G 27 - 2 G £ e 2 F 2 } , 
M n (s)\ p 2 = 0, 

M 12 ( S )| p2 = 0. (34) 

Here F n and F K are the pion and kaon decay constants respectively. 

The terms proportional to Ge are published before in fl] and, after sorting out some 
misprints there, the results fully agree. 



5.2 Next-to-leading order 

The results at next-to-leading order are a lot longer, and we decided to not include them 
all explicitely. Aq 0Q is however the full first order isospin breaking amplitude (no explicit 
photon diagrams contribute) and we therefore present it in App. [X] The expressions for 
the remaining amplitudes are available on request from the authors or can be downloaded 

M- 

The diagrams contributing are the same as in [TO], see figure there. Note however that 
vertices from C Wi can now also be from Cw2E2, £s4 also from Cs2E2 and vertices from C S2 
or £ W2 also from C E2 . 

The amplitudes at next-to-leading order can in principle include all the parameters 
from Cwa, £s4, &W2E2 and C.S2E2- The coefficients from the strong and electromagnetic 
Lagrangians are treated as known, which leaves G$,G 2 7,Nl, . . . , N± 3 , D\,D r 2 ,D\, . . . , Dj, 
-D26' • • • > -^31 and Z[, . . . , Z[ 4 as unknown. In total 41 undetermined parameters. However, 
all of these don't have to be independent in the sense that they multiply the same type of 
term, eg. m 4 ^ or e 2 m 2 . It turns out that there are 30 independent combinations, denoted 
by Ki . . . K30. See Tabled for the ones already present in the isospin limit JU] and Table |3] 
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T." 


S~1 f ATT O ATT 1 O ATT 1 ATT\ 1 /T f 1 7"Yr \ 

G 8 (iV 5 - 2A/ 7 + 2A/ 8 + N 9 ) + G 27 {~ 2 D 6 ) 


T ' 

K 2 


G 8 (N{ + N r 2 ) + G 27 (±D r 26 - J 


K 3 


G 8 {ND + G 27 [\D r 27 + lD r 28 ) 


K± 


G 27 {pi - Dl + 


K 5 


G 27 (D^o + -D31 + 2D$g) 


K G 


G 27 (81% - £% + ^0) 


K 7 


G 27 (-41% + -D29) 


Ks 


G 8 (2A£ + 4iVf + jV£ - 2iV[ - 4iVT a - 2iV{ 2 ) + G 27 (-§£>[ + fDg) 


K9 


G 8 (iV 5 r + ^ + ^) + G 27 (-| J Dg) 


Kio 


G 27 - 2L>£ - D?) 


Ku 


G27 



Table 2: The independent linear combinations of the N[ and D[ that appear in the isospin 
invariant amplitudes. 

for the new isospin breaking combinations. K 12 to Kig multiply terms including sin e, K 2 q 
to K 28 terms including e 2 and K 29 and K 30 terms including Ze 2 . 

The results proportional to Ge at next-to- leading order are also published in [TT] and, 
after sorting out some misprints there, the results fully agree. 

The corrections to the masses and decay constants including strong isospin breaking 
were also recalculated, and they were found to agree with the known results, see j201 12H] 
and references therein. 

6 Numerical Results 

6.1 Experimental data and fit 

A full isospin limit fit was made in [TU] taking into account all data published before May 
2002. One of the reasons for this investigation of isospin breaking effects is to see whether 
isospin violation can solve the discrepancies in the quadratic slope parameters found there. 
A new full fit will be done after all the electromagnetic contributions have been included 
in the amplitudes (work in progress). The data from ISTRA+ [22] and KLOE which 
appeared after JU], will then also be taken into account. 

6.2 Results with and without strong isospin breaking 

Our main result is the comparison between the amplitudes in the isospin limit and including 
first order strong isospin breaking. In Fig. |2 we show the phase space boundaries for the 
five different decays and the three curves along which we will show results for the squared 
amplitudes with and without first order strong isospin breaking. The three curves are 
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-ft 12 


^ AT r A AT r ^ AT r P. AJ r 7 AJ r P. AT r \ 
(jr 8 [6 iV 6 — 4 iV 9 — £) iv 10 — o i\ u — r iv 12 — o iv 13 j 


Kir, 

-ft 13 


(1 n /3 D r — 1 1 /1 9 dm 


J\l4 


U27 ^0 J^27 ' °/ Z - U 3l) 


^15 


n„„ ( 9/3 n r s /3 n r \ 


if 

^16 


C (O AT r _L Q AT r _L A AT r _L AT r \ 

U s [Z l\ g + 6 iV 10 + 4 l\ u + l\ l2 ) 

-a-C 1 ( Ti r i /3 Ti r \ 


ft' - 
-ft 17 


^27 


-ft 18 


Lrg (^y/ / JV 4 J + Lr 2 7 {— * L>26 ~ -^27 + ^31 J 


If 

-ft 19 


Lrg [6/Zl\ 2 ) + ^27 l-J/^-^26 + 4-^27 ~~ ^31 J 


If 

-ft 20 


(jrg \—oZj 2 + 2Z 3 — lo/O Z/4 + 0/4 Z 5 + 1/zZg + 1/Z Zy 

_i_q /J. 7 r _i_ Q /zL 7 r 4- 1 /9 7 r 4-7/in 7 r -4- 7/1 D 7 r 

-|-0/ <t Z/g to/I Z/g -f- L/ Z Zz-^q -f- / / 1U Z/-Q -f- / / 1U ^12) 




#21 


G 8 (2 2T 8 r + 2/3 ZT + 2/3 2T[ X - 2 2T[ 2 ) 


#22 


G 8 (2 Zl + 2Z r 8 + 2 Z r 9 + 4/3 Z r w + 2/3 zT[ x + 2/3 Z{ 2 + 2 Z[ 3 + 1/3 zT[ 4 ) 


#23 


^8 (-^[2) 


#24 


^8 (Zj + 3 zTg + ZJ[ + 3 Z\ x — Z{ 2 ) 


#25 


Gs (-4 Z[ - 2 ZT 2 r - 4/3 zT^ + 4 - 3/2 Z r b + zT 6 r 




-zT 7 r + 1/2 Zl + 1/2 zT 9 r + 1/3 Z r w - 1/3 Z r u - 1/3 zT[ 2 ) 


#26 


G 8 (-12/5 Zl + 3/2 Z£ - Zl + Zf - 3/2 Z r 8 




-3/2^-^-1/5^-1/5^) 


#27 


G 8 (8/5 Z\ + 4/5 zT[ x + 4/5 Z{ 2 ) 


#28 


Gs {Zl + Zl + 3/2 zT 8 r + 3/2 Z' % + Z[ + 1/2 Z\ x + 1/2 zT[ 2 ) 


#29 


Gs (3N r w + AN r n + N r 12 ) 




+G27 (D[ + 2/3D r 5 -l/3DQ 


#30 


G 8 (8/3 N r n + 8/3 iV[ 2 ) + G 27 (16/9 D r 5 + 4/9 



Table 3: The new independent linear combinations of the N£, D\ and Z\ that appear in 
the isospin breaking amplitudes. 



11 



3 




_3 I i i I i i i_ 

-1.5 -1 -0.5 0.5 1 1.5 

y 



Figure 2: The phase space boundaries for the five different decays and the curves along 
which we will compare the amplitudes. 
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5.45 


L\ 


0.38 


10- 


-3 


K x 









G27 


0.392 




1.59 


10- 


-3 


K 2 /G 8 


5.19 


10" 


-2 


Ge 


-0.4 


TT 

L z 


-2.91 


10" 


-3 


K 3 /G 8 


3.77 


10" 


-3 


5 2 - 5 


-0.576 


TT 

^4 









K4 









sin e 


1.19 • HT 2 


TT 

L 5 


1.46 


10" 


-3 


K 5 /G 27 


-4.25 


10" 


-2 


Z 


0.805 


L\ 









K§/ G27 


-1.66 


10" 


-1 


n 


0.77 GeV 


U 7 


-0.49 


10" 


-3 


K 7 / G27 


1.20 


10" 


-1 


F n 


0.0924 GeV 


LI 


1.0 


10" 


-3 


K 8 ...K n 









F K 


0.113 GeV 


K X ...K X1 









K12 ■ ■ ■ K30 










Table 4: The various input values used. 

x = 0, y = and x = y?>y. In Fig. El to Fig. [7| we then plot the different squared 
amplitudes along these curves as a function of r, where r = ±\Jy 2 + y and the sign is 
chosen in the natural way, 

f y, x = o 

r = I x/y/3, y = (35) 
I yV2, x = yV3. 

Note that for all but A^__ Q the squared amplitudes are normalized to their value at the 
center of the Dalitz plot. A comparison of the central values themselves is shown in Table |S] 
together with the decay rates integrated over the Dalitz region. In Fig.Elone can clearly see 
the thresholds introduced by the difference between and m^o induced by the breaking 
of strong isospin invariance. These thresholds correspond to a new process being allowed 
where two of the neutral pions are produced through an intermediate on shell state with 
one positive and one negative pion. 

The input values used to get the results can be seen in Table |U L\ to Lg come from 
a one- loop fit in [17j, G 8 , G27, S 2 — S and Ki to Kn from the isospin limit fit in [TU], Ge 
from and sine from ^7j. For Z we use the estimate 

z = 2F 2 e 2 ^ + ~ m ^°^ ( 36 ) 

which corresponds to the value in Table EJ As usual, \i is chosen to be 0.77 GeV. 

Very little knowledge exists of the values of K\ . . . Kn and K12 . . . K 30 , so they are set 
equal to zero (tests were also made assigning order of magnitude estimates to them). 

There are various ways to treat the masses, especially in the isospin limit case. In 
[TU] the masses used in the phase space were the physical masses occuring in the decays. 
However in the amplitudes the physical kaon mass of the kaon involved in the process 
was used and the pion mass was given by m\ = |Z)i=i,3 ?7? ' 2 r i with i = 1,2,3 being the 
three pions participating in the reaction. This allowed for the correct kinematical relation 
Si + S2 + S3 = m 2 K + Zm\ to be satisfied while having the isospin limit in the amplitude 
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Centralvalue 


Decay Rate 


Isospin limit [HI] 


Strong iso-br 


Isospin limit [TQj 


Strong iso-br 


K L -> vr u 7r vr u 
Kl — > 7r + 7r~7r° 
— > 7r + 7r~7r° 

-> 7T 7r 7T + 
fT + — > 7T + 7r + 7r~ 


6.74 • 10" 12 
7.46 • 10~ 13 


9.33- 1(T 13 
3.72 ■ 10~ 12 


6.97- 10~ 12 
7.66 • 10~ 13 


1.01 ■ 10" 12 
4.00 ■ 10~ 12 


2.65 • 10" 18 
1.63 • 10~ 18 
3.1 • 10- 21 
9.11 ■ 10~ 19 
2.97- 10~ 18 


2.74- lO" 18 
1.67- 10 -18 
3.2 • 10~ 21 
9.84- 10~ 19 
3.19 • 10" 18 



Table 5: Comparison of the central values of the amplitudes squared and the decay rates. 

but the physical masses in the phase space. The results presented used the Gell-Mann- 
Okubo (GMO) relation for the eta mass in the loops. Results with the physical eta mass 
gave small changes within the general errors given in ^UJ. In the present isospin breaking 
calculation all the pion and kaon masses were used correctly in the amplitude, but for the 
eta we again use the GMO relation, but now including isospin breaking, 

2 1 

2 / 2 1 2 2\,- L 2 lon\ 

m v = 3 ( m K+ +m K0 -m n+ ) + -m T o. (37) 

The possible lowest order contributions from the eta mass have been removed from the 
amplitudes using the corresponding next-to-leading order relation. 

Here follows a discussion of the results in somewhat more detail. In general the results 
are of a size as can be expected from this type of isospin breaking. They are of order a few, 
up to 5% in the amplitudes. The isospin breaking corrections tend to increase all decay 
rates somewhat but this can be compensated by small changes in the values of the fitted 
Ki compared to the results of [TU] . The number of significant digits quoted in Table El is 
higher than the expected precision of our results, but the trend and the general size of the 
change compared to the isospin conserving results are stable with respect to variations in 
dealing with the eta mass (physical or GMO). 

For Kl — > 7r°7r 7r° the central value of the amplitude squared increases by about 3%. 
The change in the quadratic slope is similar but the total variation over the Dalitz plot is 
small so the total decay rate increases by about 3% as well. This decay is the one which 
has most variation in the amplitude when changing how one deals with the eta mass. The 
extreme case we have found was that this effect completely cancelled the change from 
isospin violation. For this decay the amplitude as presented in the paper is also the final 
one. There are no contributions from loops with photons for this decay. It will only be 
indirectly affected when the Ki are determined from the decays involving charged particles 
which do have contributions from loops and tree level diagrams with photons. Note the 
scale in Fig. El when viewing the result. 

The squared amplitude Kl — ► 7r + 7r~7r° increases by about 2.5% with very little variation 
with the eta mass treatment. The decay rate increases by the same amount. The changes 
in the Dalitz plot slopes are similar as can be judged from Fig. 0] 

For the decay Ks — > 7r + 7r~7r° the amplitude in the center of the Dalitz plot vanishes 
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Figure 3: Comparison of Kl — > 7r 7r°7r° with and without strong isospin breaking. 
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Figure 4: Comparison of Kl — > 7r + 7r 7r° with and without strong isospin breaking. 
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Figure 5: Comparison of K$ — ► 7r + 7r~7r° with and without strong isospin breaking. 

because of the symmetries. The amplitude and the slopes increase by about 3% as can be 
seen in Fig. 

The decay K + — > 7r°7r°7r + has the largest increase. The squared amplitude in the 
center changes by about 11%. The linear slopes decrease somewhat leading to an increase 
of about 8% to the total decay rate when compared with the isospin conserved case. 

The decay K + — > 7r + 7r + 7r~ has a change of about 7.5% upwards in the center of the 
Dalitz plot and a similar change in the decay rate. The slopes decrease somewhat. 

The conclusions above do not seem to change qualitatively when we give the K\ a 
value of about 0.001 and the extra Ki a value relative to G§ and G27 of 0.01. However the 
changes induced by these values can be of the order of 10%, largest for Kl — ► 7r 7r°7r . 

It should be noted that the mentioned changes are with the values of Ki determined 
using the isospin conserving fit. A new determination including isospin breaking effects 
is planned when the diagrams with photon propagators, both at tree level and one-loop 
level, have been taken into account. At present most of the changes found can probably 
be compensated by changes in the Ki. 

7 Conclusions 

We have calculated the K — > 3tt amplitudes to next-to- leading order (p 4 ,p 2 m 2 , m 4 , p 2 e 2 ) in 
Chiral Perturbation Theory. A similar calculation was done in JU] in the isospin limit, but 
we have now included strong isospin breaking (m„ 7^ m^) and local electromagnetic isospin 
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y-o 






x = \3 y 






isospin limit 








> X 





breaking in the amplitudes . This was done partly because it is interesting in general to see 
the possible importance of isospin breaking in this process, but also to investigate whether 
isospin violation will improve the fit to experimental data. Discrepancies between data 
and the quadratic slopes from CHPT were found in ^0], and isospin breaking may be the 
cause of this. 

We have tried to estimate the effects of the breaking by comparing the squared ampli- 
tudes with and without isospin violation. The effect seems to be at a few percent level, 
and probably not quite enough to solve the dicrepancies. However, to really investigate 
this a new full fit has to be done also including the explicit photon diagrams and the new 
data (23 I2E| published after [TQj . This is work in progress and will be presented in the 
future papers Isospin Breaking in K — > 3ir Decays II and III. 
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A The amplitude for K L —> 7r u 7r u 7r u . 

We divide the M (s) function denned in Eq. (J2TJJ) as 

M (s) = M ( S )| p2+ ,|^||M 8 ^) + |^M ^( S ) + G 8 eX^) 

+ §<(*) + % M^(s) + ^ ^M 27 -( S ) + G 27 e 2 M 27 <»} • (A.l) 

The effect of G ; 8 cannot be distinguished from higher order coefficients in decays not 
involving external fields. This was known at tree level earlier and has been proven to 
one-loop in [5 . This means that terms proportinal to G' s are effectively removed. 

We have extensively used the first order isospin broken GMO relation, 

2 1 
m 2 v = -(m 2 K+ +m 2 K0 -m 2 n+ ) + -m 2 n0 , (A.2) 

in writing the amplitude in the form below. 

The explicit expressions for A, B and Bi, the finite part of the loop functions, can be 
found in many places, e.g. [2~9"j . 

The octet ones are: 

M 8 ' F (s) = 

{2L[ + 2L r 2 + Ll) (-8mJo -8/3m^„ + 8s 2 ) + (aL\ + (-8/3m>^) 
+ (N[ + N r 2 +2 N;) (m£o + 1/3 m 4 K0 - s 2 ) + N r 5 (l0/3 m 2 m 2 K0 + 2/3 m 4 K0 ) 
+N r 7 (4 m 2 m 2 K0 - 4/3 m 4 K0 ) + iV 8 r (l0/3 m 2 m^ + 4/3 m^ ) 
+N' ( - 2/3 m 2 o m 2 K0 + 2/3 m^ ) + (iV[ + 2 JVJi + iV[ 2 ) ( - 4 m 2 m 2 K o) 

+ (A(m 2 K0 ) -A(m 2 K+ )) f^^U 
\ j \fni — rrr K J 

+A(ml + ) (l/4m 2 - 1/2 m^) + A(m 2 ) (l7/24 

+A(m^ + ) (5/36 m 2 o + 25/36 777^) + A{m 2 K o) ( - 7/18 m 2 - 1/36 m^o) 
+ (A(m 2 K0 ) -A(m 2 K+ )) (5/36 m 2 + 5/18 m 2 K ) + l(mj) (l/8m| ) 
+£?(m 2 + , m 2 + , a) ( - 2 m 2 o s + mto + s 2 ) + B{m\ + , m 2 K+ , S ) (5/4 m 2 m 2 K0 
- 5/4 m 2 o s - 1/2 m 2 K0 s - 1/4 m 4 K0 + 3/4 s 2 ) 

+ 5(777 2 o,777 2 o,s) (l/2 777 2 o 777,^0 ) 

+B(m 2 r o,rn 2 K o, s) (3/8 777 2 m 2 K o — l/8m 2 K o s + 1/8 m^o) 

+B(m 2 K+ , 777^+ , S) (l/4 777 2 S + 1/4 777^0 S - 1/4 S 2 ) 
+B(m 2 K0 , m 2 s) (l/24777 2 o 777^0 + 1/8 777^0 S — 5/24777^0) 
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+B(m 2 ,m 2 , s) (l/18m^m^o) + B 1 (m 2 7T+ ,m 2 K+ , s) ( - l/2raJ m^o + 1/2 m^ ) 
+B 1 (mlo,m 2 KO , s) ( - l/4m 2 m 2 ^, + l/4m^ ) 
+B 1 (m 2 K0 ,m 2 , s) ( - 1/4 m^m 2 ^, + l/4m^- ) 

(A.3) 

M 8 ' e (,) = 

(2LJ + 2L5 + L5) (-8mt-8/3m^ + 8s 2 ) + (- 32/3 m^m^) 
+L£ ( - 40/3 m 2 m\ + 16/3 m^) + L r 7 (480 m 2 - 352 m A K ) 
+L r 8 {22Am\m 2 K - 160 m 4 K ) 

+ (iV[ + jV£ -2N; + 6N r 4 ) {-ml- l/3m 4 K + s 2 ) + JVJ ( - 26/3m 2 m\ + 10m^) 

+A^g ( - 48 m\ m\ + 40 m^) + iV 7 r (4 m 2 m^ - 4/3 m A K ) 

+N r s (22/3mlm 2 K + A/3m 4 K ) + iV 9 r (l0/3 mlm\ - 10/3 m^) 

+N{ (20 m 2 m^ - 16 m 4 K ) + (-8m 2 m 2 K ) 

+N{ 2 (68 m 2 m^ - 64 m%) + N{ 3 (96 m 2 m^ - 80 m%) 

+ o ^— 2 (^K) - ^(O) ("4) + W ( - < + 13/8 m^) 

+A(m^) (m 2 + 3/4 m^) + A(m 2 ) (5/8 m^) 

+B(ml, ml, s) (5/2m 2 m|- — 12m 2 s + 6mJ — m\ s + 6 s 2 j 

+fi(m 2 , m^, s) ( - 23/8 m 2 m^ - 7m 2 s + 5 m 4 - 3/8 m^ s + 27/8 m^ + 3 s 2 ) 

+fi(m 2 , m% s) (4/3 m 2 m 2 , - 4/3 m 4 K ) 

+B{m 2 K , m 2 K , s) ( — 2 m 2 m 2 ^ + 3 m 2 s + 4 m 2 K s — 2 m 4 K — 3 s 2 ) 
+B(m 2 K , m 2 , s) (- 11/24 m^m^ -3/8 m^ s + 23/24 m 4 K ) 
+B(m 2 ,m 2 , s) ( - 5/18 m 2 m^) +Si(m 2 ,m^, s) (5/4 m 2 m|- - 5/4 m^) 
(m 2 ,, m 2 , s) (3/4 m 2 m| - 3/4 r<) 

(A.4) 

Mt(s) = 

(3 Z 5 r - 2 Z 6 r + 2 Z£ - 3 Z r 8 - 3 Zl - 2 Z r w + 2 Z[ x + 2 Z[ 2 ) (l/9 m 2 K ) 
(-2KZ + Kl) (l/2 ml - m\) + (ff* + K r 6 ) (l/3 m 2 + 2/9 m\) 
+ {Kl + tf[ ) ( - 1/3 ml) + K{ 2 (8/3 m 2 ,) 
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4 

ml 



ml — m 2 K 



(K; - 1/2 Kl - 1/3 Kl - 1/3 Kl + 1/3 K r 9 + 1/3 K[ ) 



1 /— T/ 9\ — T/ 9 \\ / 9 ^ -, /„ 9 



? 2 (A(ml) - A{m 2 K )) ( - m\Z - 1/2 m 2 
+ o 1 yn2 %>^)(-7m^s-7/2m^s 

''* 7r "Ik ^8 

+ 4m^ + 2m^ + 3Z S 2 + 3/2^ S 2 ) 

G 8 Gs 

H ^ — 2" B ( m li m Ki s) (6ml Z s + 3ml^- s - Ami Z 

m Tv ~ m K G s 

- 2m ^-9/4Z S 2 -9/8^ 2 

Org <_r 8 

+ 2 1 2 S(m^, m^, s) (m 2 Z s + 1/2 m 2 ^ s - 3/4 Z s 2 - 3/8 §^ s 2 ) 
m%-m z K V G 8 G 8 7 

+ (A(ml) - A(ml)) (- 1/4 Z -1/8^ 
+B(ml, ml,s)(-mlZ-l/2ml^ + Zs + 1/2 ^ s 

Lr 8 ^8 

+B(ml , m 2 K , s) (7/2 m 2 Z + 7/4 m 2 ^ - 3/2 m 2 , Z - 3/4 m\ ^ 
^ , _ v .. A ,..£, S )(-l/4Za-l/8^ 

^ + l/2m^+l/4.., h 



- 3/4 Z s - 3/8 ^ s) + B(m 2 K , m 2 K , s) ( - 1/4 Zs - 1/8 ^ s) 
G 8 7 v G 8 ' 

+B 1 (ml, m 2 K ,s)(- 5/2 m 2 Z - 5/4 m 2 ^ + 1/2 Z + 1/4 m\ ^) 



(A.5) 



And the 27-plets are: 
MV(s) = 

(2L\ + 2L r 2 + LI) (8 m 4 n0 + 8/3 m 4 K0 - 8 s 2 ) + (4 + L^) (8/3 m 2 m 2 ^) 
+ (D[-3D5) (-4/3 ra>|o) +^(- H/3 m 2 o m 2 ^ - 1/3 m^ ) 

( - 1/3 m 2 m^, + 1/3 m^ ) + D r 6 (7/9 m 2 m^o - 1/3 m% ) 
+D r 7 ( - 10/3 m 2 m^ - 4/3 m 4 R0 ) 

+ (l^ 6 + 4 £>£ 7 - 3 £>£ 9 - 6 D r 30 - 6 D r 31 ) (l/3 m$ + 1/9 m 4 KQ - 1/3 s 2 ) 
;A(m 2 ,o)-A(m 2 Y+ ))(7/12mt) 

A(ml + ) ( - 5/6 m*o) + ^ * (5/3 m*o) 

IKo)(-5/6mJ„) 
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1 / 


- 

+— 


1 


m 2 „ 


- m^o 
1 


+— 

m 2 „ 


- m 2 K0 



5- B ( m l+ > m ' + > s ) ( 5 / 2 fi2 - 35 / 6 s + 10/3 mlo) 

m; - v 7 

+ o 1 2 B{ml+ , m 2 K+ , s) ( - 15/8 m 2 s 2 + 5 s - 10/3 m° ) 
— m K o 

H — 5 — - — 5— B{m 2 K+ ,m 2 K+ , s) ( — 5/8 m 2 s 2 + 5/6mt s] 

+A(ml+) (m 2 „ - 109/72 m^o) + A(m 2 n0 ) ( - 17/24 m^ ) 
+A(m^+) ( - 25/18m 2 - 25/72 m^ ) + A(m 2 K0 ) (7/18 rn^o + l/36m^„) 
+ (A(m| ) - A(m^+)) ( - 5/36 m 2 - 5/18 m^) + Z(m 2 ; ) ( - 1/8 m^) 
+B(ml+,ml+, s) (5/Gml m 2 K o - 13/6 m^ s + 2/3 m^o 

- 5/6m^ s + 3/2 s 2 ) 

+B(ml+,m 2 K+ , s) (25/6 m 2 m^ — 35/6 m 2 o s + 10/3 m^o 

- 53/24 m^o s - 1/6 m 4 R0 + 19/8 s 2 ) 
+B (mlo, mlo, s) (- l/2m^m^o) 

+B(mlo,m 2 K o, s ) { - 3/8 m^o m 2 ^, + l/8m| s - l/8m^ ) 
+B(m 2 K+ , m 2 K+ , s) ( - 13/12 m 2 s - 7/8 m 2 ^ s + 7/8 s 2 ) 
+B(m 2 K o,m 2 ,s) (- l/24m 2 r0 m^o - 1/8 m 2 K o s + 5/24 m 4 K o) 
+B(m 2 , m 2 , s) ( - 1/18 m 2 m^o) 

+5 1 (m 2 r+ ,m 2 r+ , s) ( - 19/12 m 2 m^ + 1/3 m^, - 5/12 m 4 K o) 
+ (Bi(mlo,m 2 K o, s) + 5i(m|-o,m 2 , s)) (l/4m^m|o - l/4m^ ) 

(A.6) 

M 27 ' e ( S ) = 

(2^ + 2^ + ^) (8m^ + 8/3m^-8s 2 ) 

(32/3 m 2 7/4) + L5 (40/3 m 2 m^ - 16/3 m%) 
+L r 7 ( - 480 ml m 2 K + 352 m 4 K ) + ( - 224 m 2 + 160 m%) 
+D[ ( - 76/3 ml m\ + 16 m%) + L>£ ( - 4 m 2 + 32/3 m^) 

(7/3 m 2 m\ - 9 m^) + D\ (5/3 m 2 - 5/3 m 4 K ) 
+D r 6 (73/9 m 2 m\ - 5 m^) + D r 7 ( - 22/3 m 2 m 2 K - 4/3 m^) 
+ (7 D r 2e + 28 ^ 7 - 3 D r 29 - 6 D; o - 6 Z^) (l/3 m 4 + 1/9 - 1/3 s 2 ) 

+r^\r (^«) - ( - 41/6 mj) 

m 7T _ m K 

22 



+ m2 l _ m2 B(ml, ml, s) (25/2 ml s 2 - 175/6 m\ s + 50/3 ml) 

+ „ 1 2 ff(m 2 , m^, s) ( - 45/8 m 2 s 2 + 15 s - 10 m£) 
m~ — m z K v ' 

+ 5 1 2 5(m^, m^, s) ( - 55/8 m 2 s 2 + 85/6 s - 20/3 ml) 
ml — m K v 1 

+A(ml) (41/6 ml + 43/12 m^) + A(m^) (- 41/6 m 2 - 263/24 m 2 K ) 
+A(m 2 r) )(-5/8m 2 K ) 

+B(ml, ml, s) (127/6 ml s - 38/3 m\ - 3/2 m 2 - s - 17/2 s 2 ) 
+B(ml,m 2 K ,s) (- 101/24 mlm 2 K - 59/3 m 2 s + 40/3 + 7/2 s 

- 71/24 + 61/8 s 2 ) + 5(m 2 , m 2 , s) ( - 4/3 m 2 m 2 - + 4/3 m%) 
+B(m 2 K , m 2 K , s) (26/3 ml m 2 K - 191/12 m 2 s + 20/3 

- 311/24 m^- s + 7m 4 K + 69/8 s 2 ) 

+B(m 2 K , m 2 , s) (ll/24 m 2 m|- + 3/8 m\ s - 23/24 m%) 

+B(m 2 v , m 2 , s) (5/18 ml m 2 K ) + 5 x (m 2 , m\, s) ( - 5/2 m 2 + 25/6 m^) 

+B l {m 2 K ,m%s) (- 3/4 m 2 + 3/4m^) 

(A.7) 



M 27 < e (s) = 

(2 ffj - Kl) (l/2 m 2 - m\) + (K; + Kl) ( - 1/3 m 2 - 2/9 m\) 
+ {K + X[ ) (1/3 ml) + X[ 2 ( - 8/3 m\) 

+-TZ-T ( - K l + V2 «I + 1/3 ^ 5 r + 1/3 - 1/3 ffj - 1/3 K[ ) 

+-^—1- (A(ml) - A(m 2 K )) ( - 2/3 m 2 Z) 

+ - 1 „ E(m 2 , m 2 , s) ( - 14/3 m 2 Z s + 8/3 m%Z + 2 Z s 2 ) 
ml — m A K v ' 

+ 2 - 2 £(m 2 , , s) (4 m 2 Z s - 8/3 Z - 3/2 Z s 2 ) 
+—„ T B(m 2 K , m\, s) (2/3 ml Z s - 1/2 Z s 2 ) 

+ (Z(m 2 ) - A(m^)) ( - 1/6 Z) + B{ml , ml, s) ( - 2/3 m 2 Z + 2/3 Z s) 
+fi(m 2 , rr&, s) (7/3 m 2 Z - m 2 K Z - 1/2 Z s) + S(rr&, m^, s) ( - 1/6 Z s) 
+B l (ml,m 2 K ,s) (-5/3 m 2 Z + 1/3 z) 
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-14/3m^Z(- + ln(47r) + r + l X ' 



e J 16n 2 

(A.8) 

The last term is the infinity remaining from the loops since we haven't included the 
0(G 2 7 e 2 p 2 ) Lagrangian. 
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